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Abstract A method of moment inequalities is used to derive the princi- 
ple of minimum growth rate in multiplicatively interacting stochastic pro- 
cesses (MISPs). When a value of a power-law exponent at the tail of proba- 
bility distribution function exists in a range < s < 1, a first-order moment 
diverges and an equality for a growth rate of systems breaks down. From 
the estimate of inequalities, we newly find a conditional inequality which 
determines the growth rate, and then the exponent in < s < 1. 
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1 Introduction 

Multiplicatively interacting stochastic processes(MISPs) have been investi- 
gated quite actively for past some years[T]|2]|5] | l||5] | 5] |7ll51ISini71llllll2lll3| . It is 
known that these processes generate power-law distributions at the tail of 
a probability distribution function(PDF). They are considered one of prob- 
able candidates for explaining power-law tails observed widely in natural 
phenomena [TmTB] . such as velocity distributions of inelastically colliding par- 
ticles and wealth distributions in economics. There exist two parameters char- 
acterizing the processes: a growth rate of systems and a power-law exponent. 
These are able to be calculated analytically. It is found that both parame- 
ters vary continuously depending on microscopic interaction parameters in 
a certain region. This fact shows that the power-law behavior of MISPs is 
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vastly different from that of usual critical phenomena. In critical phenom- 
ena, power-law behavior appears when the parameters are controlled only at 
a critical point. In MISPs, on the other hand, a power-law tail emerges in 
a certain region of the parameter space, which suggests ubiquity of power 
laws due to MISPs. The second difference is universality. In the former, the 
power-law exponent is independent of microscopic parameters, and therefore 
becomes an indicator of universality classes. In the latter, the exponent de- 
pends on microscopic details and universality does not hold. It is known that 
one-body random multiplicative models also show power-law tails in some 
conditions |14U16| . The dependence of the exponent on interaction parame- 
ters indicates that many-body processes is irreducible to one-body processes 
because interaction parameters are specific for many-body ones. It will be 
reported that in weighted interactions, many-body processes belong to a dif- 
ferent universality class from that of one-body ones|13|. 

The paper is organized as follows. In Sec. |2 we introduce a method of 
Fourier transform which is usually used for calculating the growth rate 7 
and the exponent s of granular gas systems 1,2,3,6 and other MISPsJSHHl 
HOIllj . In Sec. [31 we explain the principle of minimum growth rate which is 
originally found by ben-Avraham et a/.jH]. In Sec. 01 we apply a method of 
moment inequalities which is used in inelastic hard-sphere models [2] in order 
to derive the principle theoretically. In Sec El summary and discussions are 
stated. 



2 Model and method of Fourier transform 



Here, let us introduce the method of Fourier transform to derive the growth 
rate and power-law exponent flJGBK] • We consider a N— particle system de- 
fined as follows. Each particle has positive quantities Xi > [i = 1, ■ • ■ , N). 
At each time step, two particles labeled by i and j (i ^ j) are picked up ran- 
domly from N particles and interact with each other, where two quantities 
transformed into x'^x'i by the rule 

X^ — Ol.Xi I j^Xj ■ X j — j3x>i I GiX j j ^ 1 ^ 

with positive interaction parameters a, f3(> 0). When N — » 00, this process 
is described by the following master equation. 



df(z,t) 
dt 



dx dy f(x,t) f(y,t) (2) 
Jo 

x - [ 5(z — (ax + (3y)) + S(z — (fix + ay)) - S(z — x) — 6(z — y)], 

where f(z,t) is a PDF of quantities. After the Fourier transform F(k,t) = 
J °° dze lkz f(z,t), Eq. 10 becomes 

= F {ak,t)F(J3k,t) - F(k,t). (3) 



3 



In this model, the total value of the quantity is not conserved. In the long time 
limit, therefore, the PDF diverges to infinity or shrinks to zero depending on 
the value of interaction parameters a, [3. Then, the scaling relations 

£ = ze-T", ^(0 = /(«,*)e 7t , (4) 
K = fce 7 *, <P(k)=F(M), (5) 

are assumed, where 7 is a scaling parameter representing a growth rate of 
systems. Substituting Eq. (0) into Eq. J2J), we obtain 

_ $f aK \$((} K ) - #( K ). (6) 

UK 

Differentiating Eq. © with respect to k, and then substituting k = 0, we 
obtain the growth rate of systems 

7 = a + /3-l. (7) 

Performing the same calculations in Eq. J3J), we find 7 is actually the growth 
rate of the first-order moment, that is, J" °° zf(z,t)dz ~ e 7 *. Suppose <^(k) 
consists of a combination of regular and singular terms, 

W(k) = ^regular (k) + ^singular («) , (8) 



^ (*K) n 



^ n! 

n=0 



a„ + a p K p , (9) 



where p is a non-integer number in general. In this case, the small leading 
term of ^singular («) reflects the power-law tail V'CO — V? 1+S j where s is a 
power-law exponent. Substituting W(k) ~ 1 + a p K p into Eq. ©, we find a 
transcendental equation 

P7 = a p + f3 p - 1. (10) 

From this transcendental equation together with Eq. (JJJ) , we finally obtain a 
non-trivial solution p = 1) as illustrated in Fig. ^a). This solution cor- 
responds to the power-law exponent. In this way, the growth rate of systems 
7 and the power-law exponent s are calculated analytically. 



3 Principle of minimum growth rate 

In the case < s < 1, the method of Fourier transform mentioned above does 
not work because the first-order moment diverges, and thus Eq. JJJ does not 
hold anymore. Just one equation I|1U|) remains to hold and two parameters 
7 and s becomes indeterminate. ben-Avraham et aL[S] discovered that when 
< s < 1, 7 and s are determined to fulfill the following relations. 

sjmin = a s + (3 s - 1, 7 = jmin = a s In a + (3 s ln/3. (11) 

These relations indicates that the value of the growth rate 7 is equal to the 
tangential gradient of y — oP + f3 p — 1 going through the origin as shown 
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Fig. 1 Solution of the transcendental equation at s > 1 (a) and at < s < 1 (b). 
In case (b), the principle of minimum growth rate shows that the growth rate is 
given by not 7 but r y m in, and correspondingly the power-law exponent is not 5 but 
s. 



in Fig. ^b). In other words, it can be said that the growth rate 7 becomes 
minimum under the condition that Eq. (|10fl has a solution. For this reason, we 
call this mechanism the principle of minimum growth rate. To our knowledge, 
the principle has never been explained theoretically. In this paper, we derive 
the principle by applying a method of moment inequalities, instead of that 
of Fourier transform. 



4 Method of moment inequalities 

Substituting Eq. Q into Eq. ©, we obtain 

/>OC poo 

PIUp = #1 / [ Ki + K2Y - (£i) P ] . (12) 

JO JO 

where fi p is the p— th order moment of a scaled PDF ip(£,). The growth 
rate of systems Eq. {7J) is derived from Eq. H12fl at p — 1. Substituting 
the inequalities of lemma 2 in the paper by Bobylev et a/.|9], we find the 
inequalities of moments 

M a ,P,P) < < B(a,f3,p) 
7P - (dP + f3P - 1) ~ Mp ~ 7p - (a? + Pp - 1) ' 

for all p > 0, where A, B are certain bounded functions. Equation (|13fl means 
that the solution p — s of Eq. (|10|> represents the order of divergent moment 
fj, s , and thus the power-law exponent s. In this way, 7 and s are derived from 
Eqs. Q) and (|l(Jfl in the case s > 1. Therefore, this method leads to the same 
results as the method of Fourier transform did. In the case < s < 1, the 
next inequality is satisfied. 



(x + y) p < x p + yP, 



(14) 
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for all < p < 1, where x,y > 0. Substituting Eq. {H} into Eq. 1^3), we 
obtain the conditional inequality. 

P7 < a p + (3 P - 1, (15) 

for all < p < 1. It follows that when the power-law exponent exists in a 
region < s < 1, the growth rate 7 and the exponent s is determined so as 
to fulfill the next two conditions. 

1. There exist a solution p = s in < p < 1 satisfying p^ — a p + (3 P — 1. 

2. Condition pj < a p + (3 P - 1 is fulfilled for all < p < 1. 

These conditions are equivalent to Eq. Consequently, the principle of 

minimum growth rate is derived theoretically. 

5 Summary and discussions 

To sum up, we have successfully given a theoretical explanation of the prin- 
ciple of minimum growth rate in many-body stochastic processes with mul- 
tiplicative interactions. In the process of deriving the principle, we used the 
method of moment inequalities. This method made clear the correspondence 
between the solution of the transcendental equations and the divergence of 
moments as shown in Eq. (|13[) . When s > 1, the results coincide with those 
given by the method of Fourier transform. When < s < 1, Eq. 10 breaks 
down because of the divergence of the first-order moment. From the esti- 
mate of inequalities, we obtained the conditional inequality Eq. I|15fl which 
is replaced with Eq. JJ}. 

In order to obtain the principle of minimum growth rate, it is essential 
to estimate inequalities derived from the master equation of the processes. 
After performing Fourier transform, it is difficult to handle inequalities be- 
cause of the existence of imaginary parts, which is probably the main reason 
why the principle had not been explained. It can be concluded that in the 
asymptotic analysis of this kind of stochastic processes, the method of mo- 
ment inequalities are quite useful. Another application on weighted processes 
will be reported elsewhere [T^]. 
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